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We present a comprehensive theory of the circularly polarized luminescence and its dependence 
on strain in spin-polarized Si and Ge. Symmetries of wavefunctions and interactions are used to 
derive concise ratios between intensities of the right and left circularly polarized luminescence for 
each of the dominant phonon-assisted optical transitions. These ratios are then used to explain the 
circular polarization degrees of the luminescence peaks in the spectra of biaxially-strained Si and 
Ge, and of relaxed Sii-^Ge^ alloys. The spectra are numerically calculated by a combination of an 
empirical pseudopotential method, an adiabatic bond-charge model and a rigid- ion model. 
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I. INTRODUCTION 

The conservation of angular momentum during the in- 
teraction of radiation with matter allows one to study the 
angular momentum of charge carriers from the state of 
light polarization.^ ^ In semiconductors, the circular po- 
larization degree is largely set by the effect of spin-orbit 
coupling on the angular- momentum quantum numbers in 
the valence band.^ One can then use strain as a valuable 
experimental knob to regulate and validate the relation 
between the measured circular polarization degree and 
the spin polarization of carriers.^ ^ The applied strain lifts 
the energy degeneracy in the edge of the valence band, 
and it controls not only the energy spacings but also the 
mixing between hole species. Knowing the strain-induced 
hole mixing in each of the split valence bands, one can 
ascertain the spin polarization from the measured circu- 
lar polarization degree. To date, this technique has been 
widely used in direct-gap semiconductors where optical 
transitions are carried out straightforwardly.!^ 

In indirect absorption edge semiconductors, free- 
carrier optical transitions are more intricate since 
they involve both radiation-matter and electron-phonon 
interact ions .I^J^^ The photon has negligible momentum in 
comparison to charge carriers, and therefore the emis- 
sion or absorption of a phonon is needed to conserve the 
crystal momentum during electron-hole radiative recom- 
bination (between edges of the conduction and valence 
bands). Extracting information on the spin polariza- 
tion of electrons from the circular polarization degree of 
phonon-assisted optical transitions was rece ntly estab- 
lished in the case of unstrained silicon. I^E^I This inter- 
pretation allows one to quantify the spin polarization of 
electrons in experiments where strain is not a controllable 
par amet er .^^^^Ull 

In this paper we study the effect of strain on the lumi- 
nescence in spin-polarized Si and Ge. We derive concise 
ratios between intensities of the right and left circularly 
polarized luminescence for each of the optical transitions 
in these materials. These ratios depend on the phonon 
type, hole mixing, and valley position within the multi- 
valley conduction band. The derived ratios are then used 



to interpret numerically calculated spectra in biaxially- 
strained Si and Ge as well as in relaxed Sii-a^Gca^ alloys. 
Given the predictions for long spin lifetimes in strained 
Si and Ge,II3l20] ^j^^ findings of this work render optical 
transitions a viable tool in studying the spin dynamics 
when straining these materials. The dependence of spin 
relaxation on strain can be determined from the change 
in the circular polarization degree when the spin lifetime 
is comparable or longer than the electron recombination 
timescale. In this limit, the circular polarization degree is 
changed from zero to the maximal attainable theoretical 
value (that will be presented in this work). 

This paper is organized as follows. In Sec. [n| we re- 
view the theory of free-carrier optical transitions in indi- 
rect absorption edge semiconductors. After summarizing 
known spin-independent results in Ge and Si, we present 
a detailed procedure to derive spin-dependent selection 



rules. In Sec. Ill we derive expressions for the circular 
polarization degree as a function of the strain amplitude 
and type. This section also includes a brief review of 
strain effects on the energy band structure. In Sec. IV we 



confirm the analytical findings by independent numerical 
calculations that combine an empirical pseudopotential 
method, an adiabatic bond-charge model and a rigid-ion 
model. This numerical procedure is used to calculate 
the polarized spectra in biaxially-strained Si and Ge, as 
well as in relaxed Sii-a^Gca^ alloys. Section [V| is a sum- 
mary of central results. Appendix [A] includes technical 
details of the strain-dependent Luttinger-Kohn Hamil- 
tonian model. Appendix [B] deals with the moderate ef- 
fect of strain on the phonon dispersion, and Appendix [C] 
includes numerical results of optical orientation in un- 
strained Si and Ge. 



II. BACKGROUND 

The intensities of optical transitions in indirect absorp- 
tion edge semiconductors are quantified by using second- 
order perturbation theory.^^ The radiation-matter and 
electron-phonon interactions are denoted by and H^^ 
respectively, e is the light polarization vector and I is 
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the phonon mode: transverse-acoustic (TA), transverse- 
optical (TO), longitudinal- acoustic (LA) or longitudinal- 
optical (LO). The intensity of an optical transition is then 
found proportional to 
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where the angular frequency of the photon is cjq, and 
the =F sign is associated with luminescence (— ) or ab- 
sorption (+). uj£ is the angular frequency of the phonon 
emitted during the process. The initial and final states 
are given by \i) and |/), respectively. In luminescence 
and optical orientation experiments these states are typ- 
ically taken from the edges of the conduction and va- 
lence bands. The overall energy conservation implies 
that Ef — Ei = h{uJo ± oui). The sum in Eq. ^ is 
over intermediate states (|n)). The translation sym- 
metry renders crystal momentum conservation in each 
of the radiation-matter and electron-phonon interactions 
(i.e., during virtual transitions to and from intermedi- 
ate states). The interference between the first and sec- 
ond terms in Eq. ([T]) plays an important role in setting 
the amplitude of phonon-assisted optical transitions in 
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A. Spin-independent properties 

To understand the salient features of the luminescence 
or absorption spectra, one has to inspect the energy band 
structure and the phonon dispersion curves. As an ex- 
ample, in Fig. [l] we show this information for the case of 
unstrained Ge between the F and L points. The band 
structure in Fig. [ija) is from an empirical pseudopoten- 
tial method,!^ and the phonon dispersion in Fig. [ijb) is 
from an adiabatic bond-charge model.l^ The procedure 
to calculate the luminescence spectrum in Fig. [ij^c) will 
be explained in Sec.|lV[ The top edge of the valence band 
is two-band degenerate and located in the zone center 
[indicated by Fg" in Fig. [ija)]. The bottom edge of the 
nondegenerate conduction band is indicated by Lq and 
located in four equivalent zone-edge L points ([tt, tt, 7r]/a, 
[— TT, TT, tt] /a, [tt, — tt, tt] /tt and [tt, tt, — tt] / a where a is the 
lattice constant). Below we explain the origin for the 
three edge-to-edge luminescence features of the black 
curve in Fig. [TJc) as well as the origin for the difference 
in their amplitudes. 

The luminescence process is schematically described by 
combination of arrow pairs in Fig. [TJa). The vertical ar- 
rows correspond to the radiation-matter interaction and 
the arrows between L and F correspond to the interaction 
of electrons with zone-edge phonons near the L point. 
The transition from initial to final states can proceed 
via several virtual paths. From inspection, however, the 
energy proximity between the Lq and Ff states of the 
conduction band is evident (~140 meV). The phonon- 
assisted optical transitions in Ge are governed by the 
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FIG. 1: (color online) (a) Band structure of Ge along the 
r — A — L symmetry axis. The conduction (valence) edge 
state is indicated by Lq (F^). Arrows represent virtual paths 
for phonon-assisted optical transitions via pertinent interme- 
diate states. The path via F^ (L^) is shown in bold solid 
(dash) lines to represent its dominant role in setting the in- 
tensity of LA (TO) phonon-assisted optical transitions.^^ (b) 
Phonon dispersion along the F — A — L symmetry axis, (c) 
Calculated polarized luminescence due to radiative recombi- 
nation of spin-up electrons in unstrained Ge at 77 K. 



F^ intermediate state since it comes with the smallest 
'penalty' [denominator in Eq. ([T])]. Using group theory. 
Lax and Hopfield have shown that only LA modes can be 
involved in the Lq ^ F^ ^ Fg" optical transition [solid 
arrows in Fig. [ija)].^ 

In addition to the relatively strong LA spectral peak, 
the spectrum of unstrained Ge involves TO and TA spec- 
tral features.!^ When considering 'distant' intermediate 
states in Ge (the Fg , Fg and L3 ), the TO mode is also 
symmetry-allowed in edge-to-edge optical transitions,^ 
depicted by dashed arrows in Fig. flla) . However, as 
elucidated by Glembocki and Pollak,^^^the intensity of 
transition via L^ is much stronger than via F^ and F^ 
Specifically, the matrix element from electron-phonon in- 



teraction between Lo 



and Fg" 



■^5 times larger than 



that between Lq and {Fg , Fg }. As shown in Fig. [ijc) 
and supported by luminescence experiments,^^ one can 
also find a weak spectral peak from TA modes. In spite 
of the fact that edge-to-edge optical transitions with the 
TA mode are forbidden, higher order transitions between 
states near Lq and F^ do not vanish completely. The 
luminescence peak of the TA mode is visible since this 
phonon has exponentially larger population compared 
with that of other phonons. This property can be inferred 
by inspection of the L -point phonon energies [Fig. [TJb)]. 

In silicon, the absence of a single close-by interme- 
diate state renders the intensity of optical transitions 
strongly susceptible to quantum-mechanical interference. 
The six bottom edges of the multivalley and nondegener- 
ate conduction band in unstrained silicon are positioned 
~0.85x27r/a away from the zone center on the T — A — X 
symmetry axes. Invoking group theory selection rules. 
Lax and Hopfield deduced a dominant involvement of 
transverse phonons whereas longitudinal phonons suffer 
from destructive interference between valence and con- 
duction intermediate states [first and second terms in 
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Eq. ([T])]!^In experiments, however, only the TO spectral 
peak is relatively strong while the TA spectral peak is 
closer in amplitude to that of the weak LO transition!^ 
This behavior is explained by a relatively weak interac- 
tion of electrons with zone-edge TA phonons [the H£=ta 
matrix element in Eq. ([T])]. It stems from destructive 
interference between contributions of the two atoms in 
the unit cell (not to be confused with the interference 
between intermediate states). ^^^H^ 



B. Spin-dependent luminescence in Si and Ge 

To establish connection with the strain-dependent 
analysis in the next sections, we will consider the ap- 
plication of infinitesimal strain along n and calculate the 
optical selection rules for light propagation along this di- 
rection. On the one hand, this procedure provides phys- 
ical insight to the expected selection rules of strained 
crystals since it correctly classifies the zone-center eigen- 
states of the valence band. On the other hand, this 
procedure also allows us to extract the selection rules 
of unstrained crystals by neglecting the infinitesimally- 
induced energy splits in the edges of the conduction and 
valence bands. We mention that the spin-dependent opti- 
cal selection rules in unstrained Si or Ge are independent 
of the choice of light propagation (n). This indepen- 
dence is noticed from the fact that in the absence of a 
symmetry-breaking mechanism, the angular-momentum 
conservation law is invariant to the choice of the spin- 
quantization axis. 

For the general case of light propagation along n, the 
radiation operators of the left and right circular polar- 
ization vectors are 



[{px COS 6 cos (p^Py COS sin ( 
±i {py cos4>-px sin(/))] , 



- pz sin 0) 



(2) 



where Pi are vector components of the momentum oper- 
ator and {(p) are the polar (azimuthal) angles of n with 
respect to the crystallographic axes. The degenerate spin 
states of conduction electrons in the edges of the valleys 
are replaced by,'^ 



|K„^„) = cos-|K,.,^,)+sin-e*^|K,.,^,), 

|K,a„) = -sin^e-^^|K„^,)+cos^|K,.,^,), (3) 

where Kj is the wavevector in the center of the j^^ valley. 

The classification of holes states is subtle and depends 
on the form of the strain tensor. Below we discuss two 
common cases where an infinitesimal stress is applied 
parallel to the [001] and [111] crystallographic axes. The 
former (latter) represents a case where the energy de- 
generacy between conduction-band valleys in Si (Ge) is 
lifted. These strain configurations also capture the im- 
portant differences induced by diagonal and off-diagonal 



TABLE I: Representing basis functions of initial, final and 
intermediate states in the process of phonon-assisted lumi- 
nescence of spin-up electrons in germanium. Capital letters 
denote even functions, X = yz, Y = zx and Z = xy. See the 
band structure [Fig. [l]^a)] for positions of these states. 
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xyz t 


intermediate 
states 






^3" 





strain-tensor components. For example, the leading- 
order effect of [001] -strain is to couple light-hole and split- 
off states (both have Jz = ±1/2 components), while the 
heavy- hole states are left unaffected {Jz = ±3/2). On the 
other hand, of the two topmost valence bands that are 
split by [111] -strain, the one with heavier effective mass is 
described by combination of {±3/2, ±1/2} eigenstates of 
-^[111] = J • (x±y ±z)/>/3, while the band with lighter ef- 
fective mass is described by combination of{— 3/2, — 1/2} 
eigenstates of J[iii]. Detailed calculations of these well- 
established results are given in Appendix [A] 

n II [001]: The classification of the zone-center hole 
states in the presence of infinitesimal [001] -strain is sim- 
ilar to that of unstrained crystals (Appendix [A]). To ob- 
tain the intensity ratios between right and left circularly 
polarized luminescence we calculate the relative ampli- 
tudes of radiation-matter and electron-phonon matrix el- 
ements by the use of basis functions. Using the previous 
analysis and the group theory notations in Fig. [l] the 
luminescence in Ge is employed as a case study exam- 
ple. Following Sec. II A[ the conduction minima is rep- 
resented by Lg", and we take the spin- up as the 

initial state. According to this choice. Table |I] lists the 
representation basis functions of all states that take part 
in optical transitions. The notation of subscripts with 
quotation marks is to recall that the spin-orbit coupling 
is not with odd p-type orbitals but with even valence 
states. The basis functions of all possible intermediate 
states [\n) in Eq. ([T])] are listed by the direct product of 
spin- up and the basis functions of single-group irreducible 
representations. From inspection of Fig. [ija), they are 
associated with F-point conduction states and with L- 
point valence states. In Table [l| we have neglected the 
spin mixing of intermediate states caused by the rela- 
tively weak spin-orbit interaction [reflected by the small 
energy split between Fg and Fg , and between the 
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FIG. 2: (color online) Selection rules of phonon- assisted op- 
tical transitions in Si (upper panel) and Ge (lower panel). 
These rules pertain to transitions with electrons whose spins 
are oriented along the propagation direction of the detected 
light {-\-z crystallographic axis). Next to each of the indicated 
phonon modes the ratio {A+BiC+D} denotes the intensity 
ratio between right and left circularly polarized luminescence. 
Contributions from transitions with heavy (light) holes are 
denoted by A:C (B:D). 



states, see Fig. [TJa)]. However, since the summation in 
Eq. ([1]) is complete, using the concise intermediate-state 
expressions of Table [T| does not affect the final result of 
the analysis, while it simplifies the derivation. 



TABLE II: Representing basis functions of zone edge phonons. 



Mode (symmetry) 


Basis Function 


TO (L3-) 




-y)^^{^z-x-y) 


LO (L+) 
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LA (Lj) 




TA (L3+) 




-y%^^(^2z--x--y-) 



The relative amplitudes of electron-phonon matrix el- 
ements can be calculated by considering the symmetries 
of zone-edge phonons. Table [ll] lists the representation 
basis functions of these phonons. Together with the radi- 
ation operators of the two circular polarization vectors, 
^± = {Px i Wy)/V^^ the selection rules are established 
straightforwardly. For example. 
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FIG. 3: (color online) Selection rules for light propagation 
along the [111] crystallographic axis. Similar to the case in 
Fig. [2] numbers denote relative intensities of right and left 
circularly polarized luminescence of a certain phonon mode. 
These numbers should not be used to infer the relative inten- 
sities of different phonon modes. 



where Pi = {xyz\x\X) and P2 = {X \px\xy z) (dind sim- 
ilarly for their cyclic permutations) . Figure [2] summa- 
rizes the resulting selection rules in Si and Ge. For 
each of the phonon modes, we list the intensity ratios 
of the right and left circularly polarized light due to con- 
tributions of heavy-hole (HH) and light-hole (LH) final 
states. As shown in Fig. [2j light propagation along the 
z-crystallographic axis breaks the symmetry between val- 
leys in Si but not in Ge. Therefore, the selection rules in 
Si are different for transverse or longitudinal valleys with 
respect to the propagation direction, and denote 
the respective emitted circularly polarized light. 

n II [111]: We use zone-center hole states in the pres- 
ence of infinitesimal [111] -strain (Appendix A), and re- 
peat the previous analysis. The resulting selection rules 
are shown in Fig. [3j where in Ge these rules can be dif- 
ferent for the Liii valley and for the other valleys. For 
each of the phonon modes, we list the relative intensity 
ratios of the right and left circularly polarized light due to 
contributions of HH and LH states. In this strain configu- 
ration, the hole states are no longer associated with pure 
±3/2 or ±1/2 magnetic quantum numbers. Nonetheless, 
the selection rules are invariant to the choice of n when 
summing the contributions from both types of holes and 
from all valleys. This summation is relevant in unstrained 
crystals leading in the case of bulk Si to intensity ratios 
of (7+ : cr_ = 2 : 3 for either the TO or TA peaks, and 5 : 3 
for the LO peak. In unstrained bulk Ge, this summation 
leads to a 3:1 intensity ratio for either the LA or TA 
peaks, and 2 : 3 for the TO peak. 

We conclude this section with mentioning two impor- 
tant aspects. First, the 'direct-gap type' 3 : 1 intensity 
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FIG. 4: Energy shift of conduction-band valleys under biax- 
ial strain. Upper (lower) panels represent the case of Si (Ge) 
where low-energy valleys are centered on the T-A-X sym- 
metry axes (L points in the edges of the A symmetry-axes). 
Contribution to the luminescence is governed by thermal elec- 
trons that populate the low-energy valley (s). 



ratio for either the LA or TA peaks in Ge is independent 
of the valley position. The reason is that their domi- 
nant optical transitions involve electron-phonon interac- 
tion between relatively pure spin states in the lowermost 
conduction band. The second aspect relates to the en- 
ergy proximity of TO and LO phonons in Si (both are 
~60 meV with a spacing of ~4 meV). Having separate 
TO and LO spectral peaks in Si is feasible only at very 
low temperatures and high-purity silicon.^ ^ In all other 
cases the two spectral peaks merge, ^ and one should re- 
call that this single peak has a large (small) contribution 
from TO (LO) phonon-assisted optical transitions. One 
should also consider the opposite sign of their circular 
polarization degrees: —20% for TO and +25% for LO in 
unstrained Si. Due to the opposite sign contributions, 
the overall circular polarization degree of this peak is ex- 
pected to be lower than that of the less intense TA peak 
(in spite of the similar selection rules for TA and TO 
transitions). This physics elucidates the observed differ- 
ence in the measured circular polarization degree s of t he 
dominant and TA peaks in recent exp er iment s ) "'^ ^ * ^ 
Quantitative analysis of unstrained Si is found in Ref. \W . 



III. STRAIN EFFECTS 

We study the effect of strain amplitude on the selection 
rules in Si and Ge. The results of this study will be 
used in Sec.[lV|to elucidate features of the spin-polarized 
spectra that cannot be explained with infinitesimal-strain 
selection rules (Figs. |2] and |3|. Below we first describe the 
effect of strain on the energy bands, and then we derive 
selection rules in the presence of a finite biaxial strain. 
Appendix |B] includes a brief description of the strain- 
induced modification to the phonon dispersion curves. 

We start with the effect of strain on the multivalley 
conduction band. Figure [4] shows diagrams of the energy 
shifts of A and L valleys under [001] and [111] biaxial 
strain, respectively. Valleys shift up or down in energy 
depending on the angle between valley and strain axes 
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FIG. 5: Energy split and state mixing in the top of the valence 
band for two configurations of biaxial compressive strain. In 
case of tensile strain the top two bands switch positions. The 
hole state components are denoted by | J, Jn), where in the left 
(right) part n || [001] (n || [111]). The mixing between hole 
states and energy shifts are shown next to the levels of the 
strained valence bands. Details about the energy parameters 
(bcz, deo, and M^-^ and M{^^) are provided in Appendix [a| 

and on whether the strain is compressive or tensile .'^^M^ 
Evidently, the energy split between valleys suppresses the 
intervalley electron scattering. As a result, the charge 
mobility is improved. The spin lifetime is expected to 
increase more dramatic ally due to its strong dependence 
on intervalley processes .^ ^^^ 1 ^^^^ 

The effect of strain on holes is more subtle due to the 
HH and LH band-degeneracy in the F point. The trans- 
formation properties of the 4 zone-center valence states 
(including spin) belong to the irreducible representation 
of Fg". The amplitude of the spin-orbit coupling sets 
their separation from the F^ zone-center states of the 
split-off valence band. By applying strain, the Fg" energy 
degeneracy is lifted into two sets of spin-degenerate hole 
bands. Such modified bands and states can be modeled 
by a strain-dependent 6x6 Luttinger-Kohn Hamiltonian 
(Appendix [A|). The left part of Fig. [s] summarizes the 
main results for [001] -strain and the right part for [111]- 
strain. When interested in the luminescence, the most 
relevant band for optical transitions is the topmost band 
since thermal holes populate this band. When interested 
in absorption (optical orientation), other bands also be- 
come relevant if the photon energy is larger than their 
energy gap from the conduction band. 



A. Strain-dependent selection rules 

We derive spin-dependent selection rules for the polar- 
ized luminescence from strained Si and Ge. Using time 
reversal and crystal space inversion symmetries, each of 
the valence bands is two-fold spin degenerate with states 
written by 

|k,si) = K(k)|£,t)+^'£(k)|£4)]e*•^ 
|k,S2) = Yl [-&Kk)K,t)+a|(k)|€4)]e^'^-.(4) 

£=X,Y,Z 
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The states are normalized such that J2i\^i0^)\'^ + 
|6£(k)p = 1. Using [001] and [111] biaxial strain con- 
figurations, we obtain the values of and by solving 
the Luttinger-Kohn Hamiltonian. To evaluate the circu- 
lar polarization degree of dominant phonon-assisted op- 
tical transitions, we follow the straightforward analysis 
of Sec. II but we use Eq. Q to represent final states 
in the valence band. Then we find the intensities of the 
right and left circularly polarized light {I±) from which 
the circular polarization degree is extracted. 



P = 



(5) 



In the biaxial strain configuration, I± depend on the in- 
plane strain parameter 



^strain ^relax 
^^relax 



(6) 



defined by the lattice constants of the strained and re- 
laxed lattices. 



B. Silicon 

We first consider the transverse phonon-assisted opti- 
cal transitions which account for the dominant lumines- 
cence peaks in Si. We study the radiative recombina- 
tion of spin-up electrons with general hole states in the 
top of the valence band [Eq. Q]. For light propagation 
along the [001] direction, the relative intensities of left 
and right circularly polarized luminescence from x and y 
valleys read 

rJ''^^l+3{\axiY)\' + \bxiY)f)±'^^{axa*Y-bxb*Y), (7) 

where ^(x) denotes the imaginary part of x. The relative 
intensities in the ±z valleys read 



Ii^8{\az\^ + \bzn 



(8) 



The case of unstrained Si is naturally incorporated by 
assigning the values of and from the final states 
in Table [l[ Using these values, one can recover the in- 
tensity ratios between right and left circularly polarized 
luminescence in the upper panel of Fig. |2j 

[001]-biaxial compressive strain fey < 0). The low- 
ermost conduction valleys are along the x and y axes, 
and the topmost valence states are heavy holes (Figs. |4] 
and|5|. To calculate the circular polarization degree due 
to optical transitions of spin-up electrons, we substitute 
Eq. (m) into Eq. (p| and assign ax{Y) = 1(0/V^ and 
bx = oy = (Table [l|). The resulting circular polariza- 
tion degree is 



that heavy- hole states remain 'pure' in this strain config- 
uration, the — I value is compatible with the previously 
derived intensity ratio of cr+ : cr_ = 3 : 27 (see upper panel 
of Fig.|2|. 

[001]-biaxial tensile strain (e\\ > 0). The lowermost 
conduction valleys are along the z axis. The circular 
polarization degree of the luminescence vanishes in this 
case [substituting Eq. (|8| into Eq. ([5|], 



(ey > Oj — ^ 0. 



(10) 



The polarization vanishes because of the symmetry of 
electron states in the ±z conduction valleys. It is not 
caused by a valence-band related effect. 

[Ill]- biaxial strain. Our interest switches to light prop- 
agation (and spin orientation) along (x+y+z)/>/3 where 
we assign cos^ = 1/V3 and ^ = 7r/4 in Eqs. ([2|-(|3|. In 
this strain configuration, the hole states are no longer 
associated with pure ±3/2 or ±1/2 magnetic quantum 
numbers ( Appendix |A]). Using analytic al val ues of aj an d 
bi in the topmost valence bands [Eqs. ( |A20 ) and (A21)], 
we find that the circular polarization degrees are 
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(11) 
,(12) 



where Co relates to the off-diagonal component of the 
shear strain tensor. 



V3e 



Cii ± 2ci2 
cii ± 2ci2 ± 4a 



;44 



(13) 



Cij are the elastic stiffness constants which in silicon 
have values of cn = 166 GPa, Ci2 = 64 GPa and 
C44 = 79.6 GPa. Other parameters in Eq. (12) are d 



and M^^-^. The former is the shear deformation poten- 
tial of holes associated with [111] strain. Its value is 
d = —5.1 eV in silicon. Finally, 



M+i = dCo ± A; 



so 



9d^e^^ 



2o^eoAso±A|o, 



(14) 



relates to the strain-induced energy shift of the split-off 
band (see Fig. [s] for illustration). Aso =44 meV is the 
energy spacing between the topmost and split-off valence 
bands in the absence of strain. 

By assigning eo ^ in Eqs. (11) and (H^, we reach 



pr(^ii<o) 



TO.TA 



5' 



(9) 



and this value is reached if the energy split of the top- 
most valence band is larger than ksT. Due to the fact 



in agreement with the respective intensity ratios of 3 : 9 
and 5 : 3 shown in the upper panel of Fig. [3| With in- 
creasing the strain amplitude in the tensile configuration, 
the coupling with split-off states becomes stronger and 
the circular polarization degree eventually vanishes. The 
smallness of Aso in silicon is such that the drop in po- 
larization is evident already at moderate levels of biaxial 
tensile strain. Quantitatively it is seen from Eq. ( 12 ) by 
noting that M^-^-^ Mcq when dco > Aso- 

To complete the analytical study of silicon, we repeat 
the procedure for the weaker LO phonon-assisted opti- 
cal transitions. In this case, the relative intensities of 
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left and right circularly polarized luminescence are gov- 
erned by A5 intermediate states in the valence band.^ 
Using hole states from the Luttinger-Kohn Hamiltonian 
and considering the contribution from the lowermost con- 
duction valley, we get 



pr(eii<o)^o, 



Pr(e||>0) 



pr(eii>o) 

LO 



LO 



-1, 



LO 3 

~^ 5' 



- 8M+icZe<, + 32cZ2e2 



(15) 
(16) 
(17) 

(18) 



C. Germanium 

We first consider the LA phonon-assisted optical tran- 
sitions which account for the dominant luminescence 
peak in Ge. As discussed in the previous section, these 
optical transitions are mediated by relatively pure-spin 
intermediate states of symmetry. Therefore, the spin 
angular momentum is kept during the phonon-assisted 
virtual transition between the Lq and conduction- 
band states. This feature leads to similarity between 
the selection rules of direct band-gap semiconductors and 
those of LA phonon-assisted optical transitions in Ge. 

[001]- biaxial strain. This configuration does not break 
the symmetry between L valleys. We repeat the straight- 
forward analysis of Sec. |IIB| for optical transitions of 
spin-up electrons but with replacing the final states by 
Eq. Q. For light propagation along the z crystallo- 
graphic axis, we get 



LA 2^(axay - bxbl^) 



\^z\ 



(19) 



Summing the HH and LH states in Table [Tj one gets P = 
50% as expected from a 3 : 1 intensity ratio. Also similar 
to direct band-gap semiconductors is that application of 
[001] strain leads to ±100% circular polarization degree 
for optical transitions of spin-up electrons {sz = +1/2). 
It is caused by the strain-induced energy splitting of the 
valence band to heavy holes that emit only cr+ light via 
Jz = +3/2 states, or a mixture of light and split-off holes 
that emit only ct- light via Jz = —1/2 states. 

[Ill]- biaxial strain. This configuration also leads to 
±100% circular polarization degree for light propaga- 
tion along the [111] crystallographic axis. Since only 
-^[111] = ±3/2 and J[iii] = —1/2 components can take 
part in optical transitions with spin-up electrons, a com- 
plete circular polarization is guaranteed as long as these 
two components are not mixed in the topmost valence 
band. This situation is indeed valid for [111] -biaxial 
strain (Appendix \A\. The ±100% circular polarization 
degrees are compatible with the 3 : and : 1 intensity 
ratios shown in the lower panel of Fig. [3j 



We discuss the less intense TO phonon-assisted op- 
tical transitions in Ge. As mentioned before, they are 
governed by the intermediate states. For [001] strain, 
by summing the contributions from all four equivalent 
valleys, the relative light intensity reads 

/± a 4 ± 3(|az|' ± \bz\^) T 6^(axa^ - bxb^^). (20) 

In the compressive case, the topmost valence band in- 
cludes only heavy holes [Eq. ( [All )] and the resulting cir- 
cular polarization degree becomes 



pOOl 

Mr 



(^11 < 0) 



TO 



(21) 



This value is expected from the 1:7 intensity ratio of 
heavy holes shown in the lower panel of Fig. [2] The case 
of tensile strain is more involved because states of the 
topmost valence band are a mixture of LH and split-off 
states. Using Eq. (A12) and after some algebra we get 



pr(^ii>o) 



TO 



(Mo+i - 46e,)2 



(22) 



where relates to the diagonal components of the shear 
strain tensor, 



Cii ± 2ci2 



(23) 



For Ge the elastic stiffness constants are cn = 133 GPa, 
C12 = 49 GPa and C44 = 69 GPa. Other parameters in 



Eq. (22) are b and Mqq^. The former is the shear defor- 
mation potential of holes associated with [001] strain. Its 
value is b = —2.3 eV in Ge. Finally, 



^001 



bcz 



^so 



9b^el ± 26e^Aso ± A 



2 

so 



(24) 



relates to the strain-induced energy shift of the split- 
off band (see Fig. [s]). By comparing Ge and Si 
[Eqs. ( [22| ) and ([l2|], we can see that the decay to zero po- 
larization is much slower in Ge with increasing the strain 
amplitude (Aso=0.044 eV in Si and 0.29 eV in Ge). Also 
note that for application of infinitesimal strain (e^ 0), 



the ^ result by Eq. (22) is compatible to the 7:5 inten- 



sity ratio of the LH states shown in the lower panel of 
Fig-H 

The effect of [111] strain on the TO phonon-assisted 
optical transitions in Ge is slightly more involved. It 
depends on hole-state mixing and also on the lowermost 
conduction valley due to the symmetry breaking between 
L valleys. The circular polarization degrees become 



P^"(e||<0) 



Pr(6||>0) 



TO 



_9_ 

"43' 



TO 1 



(M+i - 4deo)2 



(25) 
.(26) 



By assigning eo ^ in Eqs. (25) and (|26[), we reach in 



agreement with the respective intensity ratios of 34 : 52 
and 81 : 27 shown in the lower panel of Fig. |3] 
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TABLE III: Strain-induced changes in the circular polariza- 
tion degrees of optical transitions with spin- up electrons. The 
light propagates along the strain symmetry axis. 



Phonon 



Biaxial Strain 



Polarization 



Ge LA & TA 


[001] 


compressive 


1 

2 


^ 1 




[001] 


tensile 


1 
2 


^ -1 




[111] 


compressive 


1 
2 


^ 1 




[111] 


tensile 


1 
2 


^ -1 


TO 


[001] 


compressive 


1 
5 


> 3 
4 




[001] 


tensile 


1 
5 






[111] 


compressive 


1 
5 


> 9 
43 




[111] 


tensile 


1 
5 




Si TA & TO 


[001] 


compressive 


1 
5 


, 4 
5 




[001] 


tensile 


1 
5 


-> 




[111] 


compressive 


1 
5 


, 1 

2 




[111] 


tensile 


1 
5 




LO 


[001] 


compressive 


1 

4 


^ 




[001] 


tensile 


1 

4 


^ -1 




[111] 


compressive 


1 

4 


> 3 




[111] 


tensile 


1 

4 





Table III summarizes the analytical results of this and 
the previous sections. The listed changes in the circular 
polarization degrees from unstrained to strained values, 
a ^ /3, are observable once the energy split between the 
HH and LH bands exceeds ksT (and when applicable, 
also once the energy split between conduction valleys is 
larger than ksT). In some tensile cases, indicated by 
a ^ /3 ^ for Ge and by a ^ /3 ^ for Si, the 
polarization eventually vanishes when the strain ampli- 
tude further increases. Presenting the selection rules in 
Si and Ge with ^ and respectively, is meant to recall 
that when increasing the strain amplitude the polariza- 
tion drops much faster in Si because of the proximity to 
its split-off band. 



IV. NUMERICAL CALCULATION OF THE 
SPECTRA 

To verify our symmetry-based analytical discussion in 
previous sections, we perform an independent numerical 
study of the polarized luminescence spectra. The numer- 
ical method is similar to that of unstrained Si in Ref. , 
and here we extend it to cases of unstrained Ge, biaxially- 
strained Ge, biaxially-strained Si, and relaxed Sii-a^Gca^ 
alloys. The inverse effect of optical orientation following 
absorption of a circularly polarized light is discussed in 
Appendix [C] for unstrained Si and Ge. In the next two 
paragraphes we briefly summarize the invoked numerical 
calculations. 

To obtain electronic states and energies, we em- 
ploy a spin-dependent local empirical psuedopotential 



method,^ where the strain is incorporated following the 
scheme of Rieger and Voghl^The phonon energy and dis- 
placement vectors are calculated by the adiabatic bond- 
charge model, where the strain is incorporated follow- 
ing the scheme of Eryigit and Herman. The radiation- 
matter matrix elements are calculated by the electric 
dipole approximation,^^ and the electron-phonon matrix 
elements by the rigid-ion approximation.^^ 

In the calculation of the electron-phonon matrix el- 
ements, we have interpolated the pseudopotential with 
the help of the following functional form,^ 



V{q) = aiexp(-a2(7'^) x s\n{a^q + a^). 



(27) 



To find the values of we make use of the well-known 
empirica l pseud opotential form factors at g = a/3, a/8 
and y^)2il25i4iJ Xo find a unique set of parameters, we 
also set the value of V{q = a/3/2) in Ge and of V{q = 1) 
in Si. The two values are chosen by fitting the ratio 
between intensities of different spectral peaks, Ila/Ito 
in Ge and Ito/^ta in Si, with empirical results. Since 
the dichroic behavior is generally governed by symmetry 
arguments, the selection rules are indifferent to the fitting 
procedure (see supplemental material of Ref. [9]). We 
have also verified this behavior for other interpolation 
functions. 

To make a consistent study of several strained and 
relaxed configurations, we calculate the spectra due to 
radiative recombination between spin-up electrons and 
a Fermi-Dirac distribution of holes. The electrons are 
taken from the minima of the conduction band and the 
Fermi level of the holes is positioned in the top-edge of the 
valence band. This scenario is equivalent to the lumines- 
cence that results from injection of spin-p olarize d elec- 
trons into moderately doped p-type samples.^^E^ We ex- 
pect these samples to provide optimal results in prospec- 
tive experiments. Contrary to n-type samples, the injec- 
tion of spin-polarized electrons is not masked by a large 
background of unpolarized electrons. Contrary to intrin- 
sic samples, the recombination time of minority electrons 
is greatly enhanced while their spin relaxation time is not 
severely compromised as in the case of heavily-doped p- 
type samples. 



A. Spin-Dependent Luminescence of unstrained Ge 

The numerically calculated electronic band-structure 
and phonon dispersion curves of unstrained Ge were pre- 
sented in Fig. [TJa) and (b), respectively. Figure [TJc) 
shows the spin-dependent luminescence of unstrained 
bulk Ge at 77 K. The emitted light propagates along 
the +z axis, where the red and blue solid lines denote 
the relative intensities of <j+ and cr_ polarizations, re- 
spectively. The green dashed line shows the resulting 
circular polarization degree [Eq. ([5|]. The low-energy 
spectral peak around 0.75 eV is of TO phonon- assisted 
optical transitions. The polarization of this peak is about 
—20%, in accord with the predicted unstrained value 
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FIG. 6: (color online) Calculated polarized spectra of 
biaxially-strained germanium at 20 K. The strain configura- 
tions are [001]-compressive in (a), [001]-tensile in (b), [111]- 
compressive in (c), and [lll]-tensile in (d). The red (blue) 
solid lines denote the cr+ (cr_) intensities and the green dashed 
lines are the ensuing circular polarization degrees. The anti- 
crossing of the two-topmost valence bands shown in the inset 
of (d) indicates of the strong hole-state mixing in [lll]-tensile 
strain. 



of — ^ (Table III). The dominant central peak around 



0.76 eV is of LTTphonon- assisted optical transitions and 
its polarization is positive and slightly higher than 30%. 
The polarization is lower than the predicted unstrained 
value of 50% for the LA peak, mainly because of the 
compensation from the tail of TO peak with opposite 
polarization (—20%). As will be shown below, this ef- 
fect is suppressed at low temperatures when the TO and 
LA peaks become well resolved. The spectral peak in 
Fig. [TJc) around 0.78 eV is of TA phonon- assisted op- 
tical transitions. Compared with the thermal energy, 
ksT ^ 6.7 meV at 77 K, this peak is well separated 
from the other two spectral peaks. As a result, the ef- 
fect from the opposite circular polarization degree of TO 
phonon-assisted optical transitions is diminished. This 
behavior explains why the simulated circular polariza- 
tion increases to 40% at this spectral region (close to the 
analytical value of 50%). 



B. Spin-Dependent Luminescence of Strained Ge 

We simulate four configurations of 1% biaxially- 
strained Ge. Figure |6] shows the spectra results at 20 K 
where the TO and LA spectral peaks are well resolved. 
Each of the spectra contains TO, LA and TA peaks that 
are red shifted compared with those in Fig. [TJc) due 
to the strain-induced energy reduction of the indirect 
bandgap. 

We analyze the circular polarization degrees in the 
spectra of strained Ge (dashed lines Fig. [6|. Consistent 
with the analytical picture whose results are summarized 
in Table |TTl| the dominant LA peaks show almost ±100% 
circular polarization degree. The one unique feature ap- 
pears in the [lll]-tensile configuration [Fig.pfd)], where 



the amplitude of the circular polarization degree drops 
from around —100% in the low-energy edge of the LA 
peak to nearly zero in its high-energy edge. This behav- 
ior stems from the strong hole-state mixing in this strain 
configuration as inferred from the inset of Fig. |6|d) . This 
mixing increases for states away from the F point so that 
high-energy photons can involve optical transitions with 
both types of holes, where each contributes oppositely to 
the circular polarization degree. This underlying physics 
explains the polarization drop in the high-energy edge 
tail of the LA peak in Fig. l6F 



Figure [6] also shows a weak peak in the high-energy 
part of each of the four spectra. This peak is governed 
by TA phonon-assisted optical transitions and its circular 
polarization degree is similar to that of the LA peak. As 
mentioned in Sec. II A[ both LA and TA spectral features 
are governed by conduction-band intermediate states in 
the vicinity of F^ . 

We conclude the discussion of the strained-Ge spectra 
with analysis of the TO peak. The numerically calcu- 
lated circular polarization degrees for [111] -biaxial strain 
configurations are in agreement with the analytical pic- 
ture: The polarization values in Figs. [6jc) and (d) are 
Pxo ^ —21% and Pto ^ 50% in accord with the pre- 



dicted respective results of — ^ and ^ shown in Table III 
We note here that according to the analytical picture, 
the circular polarization degree should follow the rule 
— I ^ ^ ^ wit h in creasing the amplitude of the ten- 
sile strain (Table [lll[ ). However, since the split-off band 
in Ge is relatively well separated from the top of the 
valence band, the polarization does not vanishes for 1% 
tensile strain (i.e., the — | 
is relevant). 



^ side of the selection rule 



When applying [001]-biaxial strain in Ge, the agree- 
ment between the analytical and numerical pictures is 
less successful for TO phonon-assisted optical transi- 
tions. The polarization values in Figs. [6|a) and (b) are 
Pxo ^ —57% and Pto ^ 10% wher eas the predicted 
respective results are — | and | (Table III). These differ- 
ences are attributed to the involvement of F^ intermedi- 
ate states which are no longer symmetry-forbidden. We 
can understand this behavior from the transformation 
properties of the basis function, {2z — x — y) / a/6, that 
represents the L-point TO phonons (Table [n|. Without 
strain, the matrix element of the elect ron-phonon inter- 
action between Lq and F^ states vanishes because of the 
opposite contributions of the 2z and x -\- y components. 
When applying [001] -biaxial strain, the z component is 
no longer equivalent to the x and y components, and the 
transition amplitude via F^ becomes proportional to the 
strain amplitude. This transition path has a contribution 
similar to that of the LA phonon-assisted optical transi- 
tion (50%), leading to a diminished circular polarization 
degree compared with the zeroth-order predicted value 
of —75% ( — 16.6%) in the compressive (tensile) case. 
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lated. In this configuration [Fig. [7]^d)], the circular po- 
larization degree fluctuates around a zero value because 
of the hole-state mixing in the valence band. According 
to the analytical picture, the circular polarization degree 
follows the rule — ^ ^ t ^ with increasing the strain 



amplitude (Table III ) . The proximity of the split-off band 
in Si is such that the polarization nearly vanishes for 1% 
tensile strain (i.e., the | ^ side of the selection rule 
is relevant). This behavior explains the 5% circular po- 
larization degree in the low-energy edges of the spectral 
peaks in Fig. [7|^d). 



FIG. 7: (color online) Calculated polarized spectra of 
biaxially- strained silicon at 20 K. The layout is similar to 
that of Fig. [g] The red and blue curves overlap in (b). 



C. Spin-Dependent Luminescence of Strained Si 

We repeat the analysis for the case of biaxially-strained 
Si. In the tensile configuration, we employ a stress level 
of 1% which is accessible by growing a Si layer on a SiGe 
substrate. In the compressive configuration, on the other 
hand, there are limited substrate possibilities. Here we 
employ a strain level of —0.37% which is accessible by 
growing Si on a ZnS substrate. 

Figure [7| shows the simulated spectra using the same 
layout of Fig. [6j There are two spectral peaks in each 
of the four strain configurations: a weak TA peak and a 
dominant 'TO' peak that also involves a small contribu- 
tion from LO phonon-assisted optical transitions. In the 
compressive configurations shown in Figs. [7|^a) and (c), 
the circular polarization degrees in the low-energy edge 
of the peaks are around —80% and —50%, respectively. 
These results are in agreement with the predicted val- 
ues of 



I and 

5 



- ^ shown in Table 



III 



In both configura- 
tions, the slightly larger polarization of the TA peak com- 
pared with that of the 'TO' peak stems from the opposite 
sign contributions of the nearly degenerate TO and LO 
phonon-assisted optical transitions.^ Figures [7|^a) and (c) 
also show that the circular polarization degree drops sig- 
nificantly in the high-energy edges of the spectral peaks. 
In silicon, optical transitions that correspond to these 
spectral regions involve holes that are strongly suscepti- 
ble to state mixing induced by the nearby split-off band. 

In the biaxial tensile strain configurations shown in 
Figs, [rj^b) and (d), the circular polarization degree is 
nearly zero across the entire spectral range. The zero- 
polarization is robust in the [001]-tensile configuration 
[Fig. ^h)] since the electrons populate the longitudinal 
valleys with respect to the strain symmetry axis. Con- 
sistent with the analytical picture whose results are sum- 
marized in Table [nil TA and TO phonon-assisted optical 
transitions from longitudinal valleys lead to zero polar- 
ization. When applying a [111] -biaxial tensile strain, on 
the other hand, the energy degeneracy in the conduction 
band is not lifted and all six valleys are equally popu- 



D. Spin-Polarized Luminescence from relaxed 
Sii-xGex alloys 

Depending on the mole fractions of Si and Ge, the 
conduction-band edge in Sii-a^Ge^; can be at the L point 
or around 85% alon g the A axis. The critical value of x 
is around 0.83,^^^11 which the L valleys and A valleys 
are energy degenerate. For smaller (larger) x values, the 
polarized luminescence of the alloy is expected to exhibit 
Si (Ge)-like behavior. 

Numerical simulations of the spin-dependent lumines- 
cence from relaxed Sii-a^Ge^; is based on the previous 
procedures for Si and Ge. To obtain the band structure 
and electronic states, the pseudopotential Hamiltonian 
matrices of Si and Ge are linearly combined and solved. 
Our calculation shows a critical value at Xc ~ 0.8325. To 
obtain the phonon dispersion curve and atom displace- 
ment vectors of the alloy, the dynamical matrices of Si 
and Ge are also linearly combined. 

Figure |8] shows the simulated circularly-polarized lu- 
minescence due to optical transitions of spin-up electrons 
from four compositions of Sii-a^Gca^ including the criti- 
cal composition. As before, the simulated temperature 
is 20 K and the Fermi level is positioned in the top edge 
of the valence band. The four spectra are plotted on a 
single axis to demonstrate the red shift with increasing 
the Ge content in the alloy. A breaking at the energy 
axis between 0.89 and 0.99 eV indicates that such shift- 
ing is much slower when x < Xc than it is when x > x^ 
This behavior agrees with the knowledge that when x 
decreases fro m 1 to 0, the L valleys are raised faster than 
the A valleys.123481 

We first discuss the Ge-like spectrum of Si.07Ge.93 
(lower left part of Fig. |8|. Compared with Fig. [ijc), 
the most distinctive difference is the small amplitude of 
the LA peak. On the other hand, the circular polariza- 
tion degrees of the peaks in the spectrum of Si.07Ge.93 
follow the selection rules of bulk Ge (Fig. [2|. The reason 
is that spin-dependent selection rules depend mostly on 
symmetries of the wavefunctions, whereas the amplitude 
of optical transitions is susceptible to the changes in the 
band structure. The relatively small amplitude of the 
LA peak is due to the increased energy gap between the 
Lq and valleys [larger denominator value in Eq. ([l])]. 
With reducing the Ge fraction in the alloy, this energy 
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FIG. 8: Numerical calculation of the polarized luminescence 
spectra of four Sii-xGex alloys at T=20 K. Relevant phonons 
are indicated next to the spectra peaks, with their origin of Si 
(A phonons) or Ge (L phonons). Inset: Magnification of the 
TA peak for Sio.22Geo.78- The crossing of cr+ and a- is due to 
higher order effects of optical transitions via intermediate 
states (see text). 



ure |8] shows that this luminescence is dominated by TO 
phonon-assisted optical transitions whose circular polar- 
ization degree is -20% in agreement with the derived in- 
tensity ratio of (7+ : <j_ = 2 : 3 for both Si and Ge. The 
luminescence also shows weak LA and TA peaks with 
positive circular polarization degrees of Pla ^ 30% and 
PxA ^ 17%. Pla is lower than the theoretical value of 
50% due to the overlap of the LA peak with the tail of the 
dominant but oppositely-polarized TO peak. The small 
amplitude of the LA peak, compared with pure Ge, is 
attributed to a significant larger energy gap between the 
Lq and F^ valleys. Pta is a trade-off of the two com- 
peting effects we previously mentioned for the TA peak 
in Sio.22Geo.78- 



V. SUMMARY 



gap increases due to a faster shift toward higher energies 
of the F7 conduction-band valley compared with the L 
and X valleys. 

The crossing from Ge-type luminescence to Si-type lu- 
minescence in Sii-a^Gca^ alloys takes place at relatively 
high mole fraction of Ge. The upper right part of Fig. |8] 
shows that the luminescence from Sio.42Geo.58 is already 
Si-like in nature (other than an obvious red shift in en- 
ergy). The circular polarization degrees of the TO and 
TA peaks in this simulated spectrum are around —20% 
as expected from the intensity ratio of a- = 2:3 for 
these type of optical transitions in silicon.^ However, the 
Si-like behavior is gradually modified if we keep on in- 
creasing the mole fraction of Ge toward Xc (in spite of the 
fact that electrons are still populated in A valleys). As 
an example. Fig. |8] also shows the polarized luminescence 
spectrum from Sio.22Geo.78- The inset shows the unique 
behavior of the TA peak where the <j+ and cr_ intensities 
are crossed around 1.023 eV. Here, the polarization is 
governed by the intensity ratio <j+ : (T_ = 2 : 3 in the low- 
energy edge, while it decreases to zero and eventually flips 
sign in the high-energy edge. This behavior is explained 
by the effect of intermediate states in the vicinity of F^ . 
Using a k • p expansion to describe these states, we can 
incorporate small valence-band components with ampli- 
tudes proportional to the ratio between k and the direct 
energy band-gap.^^ These components enable transitions 
with TA and TO phonons via intermediate states around 
Tj . The resulting intensity ratio of cr+ : cr_ = 3 : 1 com- 
petes with the previous 2 : 3 ratio and even overcomes 
it when the energy of the Fy valley is low enough (by 
increasing x). 

Lastly we study the luminescence of the critical al- 
loy {xc ~ 0.8325 in our simulation), where electrons 
are distributed in the degenerate L and A valleys. Fig- 



We have presented a comprehensive theory of the 
phonon-assisted optical transitions in Si, Ge and SiGe 
semiconductors. The effects of [001] and [111] biaxial 
strains on the optical properties are studied for both com- 
pressive and tensile configurations. The changes of the 
circular polarization degrees from unstrained to strained- 



crystal values are summarized in Table III The concise 
information presented in this table together with the nu- 
merically calculated spectra presented in Figs. [BJS] will be 
effective in interpretation of experiments that study the 
spin polarization of electrons from optical transitions in 
these indirect band-gap semiconductors. 

The results of our analysis also elucidate two important 
aspects that complement the analytical picture presented 



in Table HI First, in silicon the proximity between the 
split-off band and the top-edge of the valence band leads 
to significant mixing between hole states away from the 
zone center. Consequently, the circular polarization de- 
gree drops rapidly from the maximal attainable value in 
the low-energy tail of the spectral peak to nearly zero in 
its high-energy tail. In addition, the significant hole-state 
mixing leads to quenching of the circular polarization de- 
gree when applying a large [111] biaxial tensile strain in 
silicon. Second, we have identified and quantified the 
important effect of virtual transitions with the F^ inter- 
mediate states. To a large extent, changes in the circu- 
lar polarization degrees of different Sii-a^Ge^; alloys are 
caused by the fast energy shift of the F7 conduction-band 
valley when varying the alloy composition. Symmetry 
breaking by strain also leads to changes in the circular 
polarization degrees by enabling virtual transitions with 
the F7 intermediate states that are otherwise forbidden. 

This work is supported by NSF Contract No. ECCS- 
1231570 and by DTRA Contract No. HDTRAl-13-1- 
0013. 
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Appendix A: Strain-dependent Luttinger-Kohn 
Hamiltonian 

To model the hole-state mixing as well as the en- 
ergy shifts and splits in the valence band, we render a 



6x6 Luttinger-Kohn Hamiltonian.^^ When incorporat- 
ing strain, this model is also known as the Bir-Pikus 
Hamiltonian!^ The zone-center basis functions of the va- 
lence band are chosen in the order of 



X-iY^ 2Z 
-t- 



\ V2 ' VE VE ' VE VE 

Using this basis, the Hamiltonian matrix reads 



X-iY 



V2 



i, 



X+iY 



V3 V3 



X-iY 



/ P + Q 


-S 


R 





s 


V2R 


\ 


-S* 


P-Q 





R 


-V2Q 






R* 





P-Q 


s ^ 


lis* 


V2Q 







R* 


S* 


P + Q - 


V2R* 


S" 




S" 


-V2Q 




-V2R P 


+ Aso 







y V2R* 




\/2g 


s 

V2 





P + Aso 


J 



(Al) 



(A2) 



where Aso is the split-off energy and the other param- 
eters include wavevector components {ki) and strain- 
tensor components (e^j). 



2mo 
2mo 



^kl- 2kl) - -{e^^ + Cyy - 2e^^), 



H ^^(^xx ^yy^ ^d^xyi 
S = ^^^^^^^3(^0; - iky)k^ - d{e^^ - icy^). (A3) 

{71772,73} are the Luttinger parameters, and ay is the 
volume deformation potential giving rise to a uniform 
strain-induced energy shift of the valence bands, b and 
d are the shear deformation potentials giving rise to the 
energy split in the top of the valence band when applying 
[001] and [111] strain, respectively. The analysis of the 
circular polarization degree, discussed in the main text, 
is almost independent of the values of these well-known 
empirical parameters. In the infinitesimal strain picture 
{cij 0), the only useful information for deriving the 
selection rules is that the shear deformation potentials 
{b and d) are negative in Si and Ge. In the finite strain 
picture, we also consider the proximity of the split-off va- 
lence band in Si (Aso = 44 meV). This proximity leads to 
a measurable variation in the circular polarization degree 
with a moderate increase of the photon energy. 



Following the discussion in the main text, the solution 
of the Hamiltonian in the zone-center {k=0) is sufficient 
in order to derive the strain-dependent optical selection 
rules. Below, we present the results for the biaxial strain 
configurations studied in this paper. 



[001] biaxial strain 

This configuration is found when thin layers are grown 
on lattice-mismatched [001] substrates. The strain tensor 
is diagonal with components. 



^yy 



2Ci2 
Cll 



0, 



(A4) 



where ey was defined in Eq. (|6|, and Cij are the elastic 
stiffness constants with values provided after Eq. (13) for 
Si and after Eq. (I23|) for Ge. The solution of Eq. ( A2|) 



for zone-center states is simple and can be compactly 
expressed by defining the constants. 



001 



+ Aso ± A/962e2 +26e, Aso + A|o , (A5) 



where = — \ {^xx ^^yy)- Using these constants, the 
zone-center energies are 



El — ay{e^ 

E2 = CLvi^x 



.) - be. 



- ^yy -\- ^zz) -\- 2^*'^001' 



E3 — Civi^xx + ^yy ^zz) + 



001- 



(A6) 
(A7) 

(A8) 
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Recalling that b is negative in both Si and Ge, the energy 
of the topmost band is Ei in biaxial compressive strain 
(ey < 0) and E2 in biaxial tensile strain (ey > 0). Apart 
from the mutual energy shift, the left-hand side of Fig. |5] 
shows the relative arrangement of £^1, E2 and Es in the 
compressive configuration. 

The circular polarization degree of the luminescence is 
governed by optical transitions with holes of the topmost 
band (ground state). In the compressive configuration, 
the topmost zone-center eigenvectors are from heavy-hole 
states, 



Vg7(e|| <0) = [1,0,0,0,0,0]^, 

VgHH^w <o) = [0,0,0,1,0,0]^. 



(A9) 



In the tensile configuration, the topmost zone-center 
eigenvectors are a mixture of light holes and split-off 
holes, 



nT(^ii>o) = 



1 

N 



nT(^ii>o) = ^ 



0,0,^,0,0,-6e, 
0,^,0,0,-6e„0 



(AlO) 



N=y |Mqq^ + 6^e^ is a normalization factor. The com- 
ponent of the split-off holes is proportional to the strain 
amplitude. In the case of infinitesimal strain, we assign 
Mqq^/N = ^/S and = 0, and get that the eigenvec- 
tors are from pure light-hole states. Using the expansion 
in Eq. Q, we find the coefficients from which we derive 
the selection rules. For example, the coefficients of the 
spin-up state are 



ax 



1 i 

"71' '^""71' 

az = bx = by = bz = 0, 



(All) 



in the [001]-compressive configuration, and 

1 Mq+i _ J_be^ 
4v^ iV 73 AT ' 

4V3 N VSN' 
1 1 be. 



ax 
ay 

bz 

az 



2Vi N 73 AT 
bx=by = 0, 



(A12) 



in the [001]-tensile configuration. 



[Ill] biaxial strain 



of the strain tensor are 



^yy 



4C44 



Cii + 2ci2 + 4C44 " 

Cii + 2ci2 

Cii + 2ci2 + 4C44 



e||. (A13) 



The solution of Eq. (A2) for zone-center states can be 



compactly expressed by defining the constants 

M±i = deo+Aso±A/M2e2+2o^eoAso+A|o, (A14) 



where Cq = {exy + e 



yz 



c)/a/3- Using these constants. 



the zone-center energies are 

E2 = Civ{exx + ^yy ^zz) + 
E3 = CLv{^x 



1 



(A15) 
Mill, (A16) 



1 



^yy 



e..)+^eo--M+i. (A17) 



Recalling that d is negative in both Si and Ge, the energy 
of the topmost band is Ei in biaxial compressive strain 
{cq < 0) and E2 in biaxial tensile strain (e^ > 0). Apart 
from the mutual energy shift, the right-hand side of Fig. [5] 
shows the relative arrangement of Ei^ E2 and Es in the 
compressive configuration. 

The order of the two-topmost bands is set by the sign 
of ey where we recall that d is negative in both Si and 
Ge. Apart from the mutual energy shift, SayCd, this in- 
formation is shown in the right-hand side of Fig. |5] for 
the compressive configuration (ey < 0). 

The circular polarization degree of the luminescence is 
governed by optical transitions with holes of the topmost 
band (ground state). In the compressive configuration, 
the topmost zone-center eigenvectors are a mixture of 
heavy and light holes. 



nT(^ii<o) = 



n!"(eii<o) 



1 i 

71 '71'' 



1 



V6 



,0,0 



(A18) 



In the tensile configuration, the mixture is between all 
species of holes. 



nT(^ii>o) 



^;,"(^ii>o) = 



1 

TV 



2V6 



4V2 



M+1,0, 



4V6 



M+i,0,deo 



1 

TV 



1 



4^6 



4V2 



2^6 



(A19) 



This configuration is found when thin layers are grown 
on lattice- mismatched [111] substrates. The components 



where the normalization factor is A/" = y l^ni + <i^e^. 
Due to the off-diagonal nature of the strain-tensor one 
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cannot associate the ground states with pure ±3/2 or 
zbl/2 magnetic quantum numbers. It is true even when a 
vanishingly smah strain is apphed [assign My^^/N = ^/S 
and 60 = 0] . 

In the infinitesimal strain picture we can still classify 
the topmost band by its effective mass: HH in the [111]- 
compressive configuration and LH in the [lll]-tensile 
configuration. To better understand the spin-dependent 
selection rules of these bands, it is convenient to make 
use of the angular- momentum eigenstates of J[iii] = J-n, 
where J is the vector of spin-3/2 angular momentum ma- 
trices and n = (x + y + z)a/3. From the four eigenstates, 
{|=b|)iii,|±|)iii},we get the following information when 
a vanishingly small strain is applied. The composition of 
the spin-degenerate ground states is from J[iii] = +| 
and J[iii] — +^ components in the [111] -compressive 
configuration (heavy holes), and from J[iii] = — § and 
^[111] — ~\ components in the [lll]-tensile configura- 
tion (light holes). We make use of this information when 
deriving the selection rules for optical transitions with 
spin-up electrons (5'[iii] = ^ in this configuration). 

Finally, using the expansion in Eq. Q, we find the 
coefficients from which we derive the selection rules. For 
example, the coefficients of the spin-up state are 



1 



-1 



ax 
bx 



2^3' 
1 + i 
2^3' 



ay 

by 



2^/3 ' 



71' 



az 

hz = 0, (A20) 



for [111] -compressive strain, and 



ax 



ay = 



4V3 N 

4^3 N ' 
for [lll]-tensile strain. 



1 dco 

i dco 
71 A^' 



1 M+i 



az 



bx = 
bY = - 
bz = - 



1 <ii 
4\/3 ■ 
1 



4\/3 N 
1 dtp 
\/3 N ' 



(A21) 



Appendix B: The effect of strain on the phonon 
dispersion 

As explained in the main text, changes to the mag- 
nitude, position and polarizations of the spectral peaks 
are mainly caused by the strain-induced degeneracy lift- 
ing and hole-state mixing.' ^^^^^ ' Changes in the phonon 
dispersion have a lesser effect. In the strained crys- 
tal, th e bon d length and bond angle between atoms 
change .'^^^ However, as long as the strain levels are not 
very large, one can still assume an harmonic oscillation 
of atoms about their new equilibrium position. We have 
generalized the adiabatic bond-charge model to account 
for strain. Figure [9] shows the phonon dispersion of 
biaxially-strained Ge in the [111] compressive configura- 
tion. From inspection of the lifted energy degeneracies. 



one can infer that the resulting broadening of the TA 
and TO spectral peaks in the luminescence of strained 
crystals will be negligible. Similarly, the expected energy 
shift of all spectral peaks is minute. 




FIG. 9: Phonon dispersion along inequivalent A symmetry 
axes in Ge when applying [111] biaxial compressive strain of 
1%. The two- fold degeneracy of transverse modes is lifted if 
the phonon wavevector and strain directions are nonparallel. 



Appendix C: Optical Orientation of Si and Ge 

Optical orientation is the inverse process of circularly- 
polarized luminescence. It involves the transfer of angu- 
lar momentum from a circularly polarized light an elec- 
tronic spin system. Here we focus on the phonon-assisted 
optical transitions across the indirect band-gap of un- 
strained Si and Ge. We will not cover the absorption in 
the direct-gap spectral range which for the case of Ge is 
very effective in generating spin polarized electrons (due 
to the energy proximity between the Vj and L conduction 
valleys). I^SMl 

The spin polarization of photoexcited electrons is stud- 
ied following absorption of a circularly-polarized light 
that propagates along the +z axis. To interpret the sim- 
ulated spin-polarization, we make use of Eq. ([T]) and the 

'2] Using similar 
we calculate the 



spin-dependent selection rules in Fig. 
numerical procedure to that of Sec. |IV 
spin-resolved absorption coefficient a following light ex- 
citation with (7+ polarization, a depends on the density 
of states in the conduction and valence bands. In bulk 
materials it shows quadratic dependence on the photon 
energy. In Fig. 
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we show the calculated value of a^/^ 
versus the photon energy for Si and Ge. 

At low temperatures, the free-carrier absorption edge 
of Ge amounts to the sum of the indirect band-gap energy 
and the energy of the L -point TA phonon [Fig. [Toj^a)]. 
Using the selection rule of this transition, the ratio be- 
tween the excited spin-up and spin-down electrons should 
be 1:3, leading to spin polarization around —50%. The 
onset of LA phonon-assisted optical transitions leads to 
enhanced absorption shown by the larger-slope region in 
Fig.jlOja). Since the LA and TA selection rules are iden- 
tical, the injected spin polarization of electrons does not 
change appreciably. The onset of TO phonon-assisted 
optical transitions at even higher photon energies leads 
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FIG. 10: Numerical calculation of a^^^ for spin-up and spin- 
down electrons by absorbing cr+ light in Ge (a) and Si (b). 
Note that in Si we resolve the contributions for a^^^ from 
transverse and longitudinal valleys where the curves are ver- 
tically shifted for better illustration. The dashed green lines 
denote the spin polarization of electrons from all valleys. 



to smaller spin polarization because of the opposite se- 
lection rule (3 : 2 versus 1:3). Finally, when the photon 
energy approaches the direct-gap energy (but below its 
absorption edge), the LA phonon-assisted optical tran- 
sitions intensify and the spin polarization of electrons 
increases. 



Figure 10 'b) shows the spin-resolved and valley- 
resolved value of Oi^l^ for Si. The absorption edge 
amounts to the sum of the indirect band-gap energy and 
the energy of TA phonons along the A-axis near the zone- 
edge X point. The spin polarization of electrons from 
all valleys is around 20% at the absorption edge. The 
onset of the dominant TO phonon-assisted optical tran- 
sitions leads to enhanced absorption. Most importantly, 
the spectral window for efficient optical spin injection in 
Si is of few tens of meV because of the proximity of the 
split-off band. 

According to the selection rule in Fig. |2j the ratio 
between spin-up and spin-down electrons excited to x 
or y valleys in Si should be 2 : 1 for transverse phonon- 
assisted optical transitions. This information is shown 
in the low part of Fig. 10 'b) where the blue solid line 



(spin-up) is always higher than the red solid line (spin- 
down). For electrons excited to the z valleys, the spin 
polarization can only come from the relatively weak LO 
phonon-assisted optical transitions for which the ratio is 
1:3. In calculating the total spin-polarization value we 
have assumed that all valleys are equally populated. The 
spin-momentum correlations and momentum alignment 
during excitation were neglected.^^^ These effects may 
introduce quantitative but not qualitative changes (due 
to valley repopulation). 
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